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We investigate the Dvali-Gabadadze-Porrati (DGP) model corrected by higher curvature brane 
terms. We show that these corrections have a dramatic impact on the spectrum of the model at 
the linearized level. Owing to the presence of higher derivatives in the field equations very massive 
ghost excitations with mass of order of Planck mass are generated in the ordinary branch of the 
model. These excitations describe an instability of Minkowski vacuum with time-scale of order of 
the Ifubble time H^^ . At large distances these tachyonic excitations are expected to decouple from 
brane-locaUzed matter. Our modified DGP model represents therefore a very promising framework 
for solving of the cosmological constant problem, in which Planck-scale physics is responsible for 
the elementary excitations driving the accelerated expansion of the universe, but the time-scale of 
the instability is settled by gravitational physics at large scales. 



Gravity models that allow for large-distance deviations from standard General Relativity (GR), the so-called Mod- 
ified theories of gravity [H, 0, 13, 0] > have been widely investigated in recent years. The main motivation behind this 
interest is the hope that the accelerated expansion of the universe [1, 0, 0, 0| could be explained and the cosmological 
constant problem solved by very large distance, of order r 3> Tc = H^^ ^ modifications of General Relativity without 
postulating nonbaryonic forms of matter such as dark energy. From a purely theoretical point of view Modified 
Gravity represents a way to circumvent Weinberg's no-go theorem on the "old" cosmological constant problem 0, 
and to shift it from the realm of short distance particle physics to that of large distance gravitational physics. More- 
over, from the phcnomenological point of view there is plenty of room for such infrared modifications of gravitational 
physics. Our experimental knowledge of gravity is in agreement with GR but is limited to distances between say 
lO^^cTO- lOOMPc. 

One of the most promising scenarios for modifications of GR in a rclativistic and general covariant framework is the 
Dvali-Gabadadze-Porrati (DGP) model % The DGP model works in the context of a brane- world scenario (see Ref. 
[ill] for a early proposal), in which bulk non-local effects modify the large distance behavior of gravity on the brane. 
In particular, the DGP model allows in the so-called self- accelerating branch for tachyonic excitations describing an 
instability of Minkowski vacuum with a time-scale of order t/ = Tc. Such an instability could represent a solution to 
the cosmological constant problem. Unfortunately, it was soon realized that the tachyonic excitation is a light ghost 



of mass m ~ tod — r'^ gravitationally coupled to brane- localized matter |12l.ll3j. From particle physics point of view 
the presence of such a light ghost is simply disastrous. 

There is no compelling reason for having t/ ~ . We would rather expect the cosmological constant problem to 
be solved by a conspiracy between short-distances and large-distances physics. Motivated by these considerations, in 
this letter we explore the possibility of a brane- world scenario in which the the time-scale of the tachyonic instability is 
still of order rc ( i.e it is determined by large-distance gravitational physics), but its mass is of order of 4D Planck-mass 
Mp. Introducing terms quadratic in the curvatures in the brane part of the DGP action we show that at the linearized 
level the model allows, in the ordinary branchy for tachyonic excitations describing an instability of Minkowski vacuum 
with timescale rc- Although these excitations are ghosts, they are heavy with mass ~ Mp and at large distances 
they are expected to decoupled from brane-localized matter. 

In the DGP model all known interactions, except gravity, are confined on a zero-tension (1 -I- 3)-brane that is 
embedded in a five-dimensional (5D) infinite- volume space-time (the bulk). Although 5D brane-induced gravity 
cannot solve the old cosmological constant problem [l^ . it can be used to deepen our understanding about modified 
gravity theories. The DGP action is 

S^M^ j d^xVC (^'i? -I- AfJ, J d^x^f^R + bound, terms , (1) 

where M, X, G and '^'i? are bulk quantities, while Mp, x, g and R are referred to the brane. In Eq. ([T]) boundary, 
Gibbons-Hawking, terms should be taken in account to warrant the correct equations in the bulk. Matter fields are 
also localized on the brane and they are omitted here for simplicity. 

The relevant new ingredient in Eq. ([1]) is the Einstein-Hilbert term on the brane. This 4D kinetic term, as well as all 

the possible Lorentz invariant terms, can be induced via loop corrections by brane-localized matter. 

Using the notation tjab = (H ), with capital indices running on 0,1,2,3,5 and Greek indices running on 
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0, 1, 2, 3, the field equations stemming from the action ([T]) are, 

Gfs = ^'5(brane)<5^^<5^(r^. - 2M|.G^.), (2) 

where is the Einstein tensor built from the 5D metric, g\g^ and G^v is the 4D Einstein tensor built from the 
induced 4D metric, = d^X^dvX^ g^^^. X^{x'^) defines the brane position in the bulk. We consider a static 
brane, located at y =: = 0, endowed with a Z2-symmctry. Eqs. ([2]) are extremely difficult to solve, only symmetric 
solutions are known. Deffayet found the modified Friedman- Robertson- Walker cosmological solutions [l^l , while some 
approximation for the Schwarzschild solution has been also derived 3, 17| . 



Relevant physical information about the model (spectrum of excitations, stability) can be gained at the linearized 
level by expanding the metric around the Minkowski vacuum: = tjab + hAB- The {fiv} component of Eq. 12]) in 
the harmonic gauge, d^hAB — ^^sh^, becomes 

[M^dcd'' + Ml6{v)d^d'^] h^, = -Siy) {T^, - Ij^^.T) + 5{y)Mld^,d,h^^ . (3) 



The structure of the propagator in the Fourier space is 
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V = o)r^''(p) = — , (4) 

Mp p'^±mDP 

where p = yjp^ is the square root of the Euclidean 4-momentum, and mjj — 2M^ /Mp is the DGP crossover parameter. 
The different signs in Eq. ^ are related to different boundary conditions at y — > oo and to the brane embedding 
in the 5D bulk. They describe different branches of solutions. The sign refers to the ordinary branch and the 
"— " sign to the self- accelerating branch. The spectrum of linear excitations near the Minkowski vacuum can be read 
out from the pole structure of the propagator (j4]) . In the ordinary branch we have a massless excitation, which does 
not propagate physical degrees of freedom (corresponding to a pole with zero residue) and a spin 2 resonance of mass 
771 — niD- In the self- accelerating branch we always have the massless excitation and a tachyon of mass to_d. 

The tensorial part of Eq. Q is the same as that of Pauli-Fierz (PF) theory for a massive graviton. In fact, from 
a 4D point of view, gravity in the DGP model is mediated by a continuum of massive Kaluza-Klein modes, with 
no normalizable massless graviton in the spectrum. This is a consequence of the presence of a flat and infinite size 
extra-dimension. We have therefore a 5D (or a massive 4D) graviton propagating 5 on-shell degrees of freedom: a 
helicity-2 state which is the analog of the GR massless graviton, a helicity-1 state which does not contribute to the 
propagator because of the conservation of the stress tensor and moreover there is a helicity-0 state. The presence of 
the latter state is the cause of the van Dam-Veltman-Zakharov (vDVZ) discontinuity [l^l- This is related to the fact 
that in the limit ?7i£) — > of the linearized approach, Newton potential is not exactly recovered but the Newtonian 
constant Gn has to be renormalized. 

The static potential of the DGP model behaves therefore like V{r) ~ 1/r for r ^ rc — 1/mo, while it behaves 
as a truly five dimensional Newton potential, Vir) ^ l/r^, for r ^ rc- It is usually required rc ~ in order to 
allow for IR modifications of GR. 

Although the vDVZ discontinuity makes the linear approximation of the DGP model problematic, it has been 
argued that nonlinearities may play a fundamental work, restoring solutions sufficiently close to those of GR at solar 
system scales [H, HO, IH . 

From a phenomenological point of view the self-accelerating branch is the most interesting one, because it allows 
for cosmological de Sitter solutions even without a 4D cosmological constant J^]. In fact the tachyon corresponds 
to an instability of the Minkowski vacuum with an exponential growth time, e™"*. This is a very appealing feature, 
since the background is modified with a curvature term of order ~ rn^. The new parameter, mo, has to be fine-tuned 
to describe the present acceleration of the universe, thus it may seem that there is no advantage in using the DGP 
model instead of introducing a cosmological constant. But, in the DGP case the acceleration of the universe is a truly 
dynamical effect, which is a consequence of the modified field equations. 

Unfortunately, the DGP model has been found to suffer from different problems. There is a strong coupling problem 
emerging when one tries to extend the results beyond the linear approximation l^, 23 1. At the first non linear order 



the propagator gets divergent contributions. But the most dangerous issue is the presence of ghosts in the spectrum. 



It has been shown |12l.ll3l| that the tachyon in the self-accelerating branch is a ghost of mass mo, i-e. it has a negative 
kinetic energy term. The presence of a light gravitational ghost excitation that couples to brane-localized matter 
makes the 4D description meaningless and seems to rule out at least the self-accelerating branch of the DGP model 
as explanation for cosmic acceleration [ist . 
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On the other hand the ordinary branch is ghost-free, but in this case the Minkowski vacuum is stable and no de 
Sitter solutions exist without introducing a positive cosmological constant on the brane. This is however strictly true 
only if we consider a brane action of the Einstein-Hilbert type. The situation could change if we allow for other 
general covariant terms in the brane action, e.g terms depending non-linearly on the 4D curvature tensors. Motivated 
by this arguments, we consider a DGP model with also induced second order curvature terms on the brane (28j : 

S = [ (fX^/GR(r^) + / d^xy/^ {mIr + aR^.uR'"' - PR^) , (5) 



H5) 

where a, (3 are dimensionless real parameters. 

The DGP model is recovered ioi a — P = 0. It is worth to notice that, due to the 4D Gauss-Bonnet relation the 
quadratic terms in ([5]) are the most general second order terms that can be added on a (l-|-3)-brane. There are several 
motivations for considering these corrections to the DGP model. First, the same quantum corrections which induce 
the Einstein-Hilbert term on the brane, can also induce higher order terms via higher loop corrections. Second, the 
introduction of these contributions could shed new light on some pathological aspects of the DGP model. Third, the 
action ([5]) represents the most general contribution at linearized level of higher curvature brane terms. This is because 
third (or higher) order curvature terms do not contribute to the linearized field equations (their contribution to the 
field equations are of 0{hf^j,)). Last but not least, second order curvature terms in the action ([5]) are suppressed by 
Mp^ powers of 4D Planck mass. They represent short-distance brane corrections to the DGP model. The action ([5]) 
represents therefore a nice framework to investigate simultaneously IR and UV modifications of GR. This is the most 
promising setting for solving the cosmological constant problem. 

For M = the theory (O reduces to the higher derivatives gravity theory first proposed by Stelle [1^ . This model 
has some interesting features: it is renormalizable and it contains massive spin 2 modes. Nevertheless, it is well known 
that in general such higher derivatives theories are pathological, due to the unavoidable presence of ghosts. In Stelle's 
model the ghost has a mass ^ Mp/^/a and can be therefore exited at high energy scale. This makes the theory 
inconsistent at least if one requires unitarity. 

In the framework of brane-world models the unitarity requirement can be reformulated in a different way. Actions 
like 111 or ([5|) have to be considered as low-energy effective descriptions of some fundamental unitary quantum theory 
of gravity, such as string theory. At short distance, say of order of Planck length, actions H]) and (O are not 
longer valid and the infinite series of higher order terms should be taken into account. For instance, in the complete 
fundamental theory higher curvature terms can contribute in some form protected by topological invariance (as in 
the Gauss-Bonnet case) such that unitarity is restored. 

Thus, as far as we are interested in the large distance, IR behavior of gravity described by the model (O, we can 
neglect the presence of ghosts as long as they have a mass ~ Mp which is above the UV cutoff of our effective theory 
and they do not couple to brane-localized matter. 

The modified DGP field equations derived from the action ([5|) are, 



with 

V = {a-2l3)V^,V,R-aV^V[}R,,,+ 

- (|-2/3)5^,V'^V0i?-f 2ai?"X-/3 + 2/3i?i?^. + i<?M-(ai?^''i?p.-/3i?') (7) 
Following the original computation of 1], the linearized equations for the {fiv} component read 

[M^dcd'' + Ml5{y)d^d'' + adpdPd^d'^] h^, = -6{y){T^, - ^v^^T) + S{y) [M].d^d, + ad^d^d^.d''] /igs- (8) 

Notice that the linearized field equations ^ do not depend on the parameter /3. This is due to the fact that in the 
harmonic gauge the 4D linearized Ricci scalar is identically zero. This means that apart from the R term, only the 
(-R^i/)^ term contributes to the propagator. In particular the introduction of /(i?) terms on the brane does not change 
the DGP propagator (gl). 
For a, /3 the propagator turns out to be 

1 / T T^^'^ — -T^ \ 



4 



The tensorial part of the propagator is identical to that of the DGP model and to that of massive gravity. But in our 
case the excitation spectrum of the theory is much richer and it depends not only on the two branches of the theory, 
but also on the parameter space {mjj, a, Mp). Moreover, owing the the higher derivative terms in Eq. ([5]) we expect 
the presence of ghosts in the spectrum. 

As usual, the spectrum of linear excitations can be inferred from the pole structure of the propagator Q . In the 
[Re(p), Im(p)] plane the propagator © has in general 4 complex poles. One is located at = 0. As in the DGP model 
the residue of this pole is zero, consistently with the absence of a normalizable 4D massless state in the spectrum. 
For a ^ the other three poles can be found solving the algebraic third order equation 



Ml 



P 



a 



-P± 



mnMp 



0. 



(10) 



Defining ag = (mf") ^^'^ different regions in the parameter space: (i) for a > and a < —ao there is one 
real solution and two complex conjugate solutions, (ii) for — ao < a < there are three real solutions. Notice that 
the case a < — ao is phenomenologically highly suppressed since ao 3> 1. In any realistic situation we are left with 
a > and a small a < 0. These two regions are discussed in detail below. 

In the following we will focus on the ordinary branch of the theory corresponding to the + sign in equation ([9|). 
Note that, due to the absence of terms in Eq. pO]) the sum xi + X2 + X3 of the poles is identically zero. 
For a > the solutions of Eq. (fTO|) are 



mi 

xi ^ -mi X2,3 = — ±m2i , 



where we have defined the "effective" masses 



mi 



m2 



M 



a3 




a/3M 

2a^ 



l + x 



l + x 



(11) 

(12) 
(13) 



and X — —. In any phenomenologically viable model we have naturally a ~ 0{\) and x is an extremely small quantity 



{x ^ 10 with reasonable values Mp 
the masses become 



lO^^GeV, viD ~ i?o ~ 10 ^^GeV and a — 1). With this approximation 



mi 



m2 



rriD 




(14) 
(15) 



As expected there are two independent energy scales in the model: m^j , which is the scale arising in the DGP model, 
and 171$ = ^^-j which is exactly the mass of the massive 4D graviton found by Stelle in the framework of gravity 
theories with higher derivatives. The pole structure of the propagator is shown in Fig. ([1]). 

Interesting novel features arises from the investigation of pole structure represented in Fig. ([1]). Similarly to the 
ordinary branch of the DGP model, there is a resonance state, corresponding to the pole p = —mi ~ —rno- This 
particle couples to brane-localized matter and does not lead to any instability of the Minkowski vacuum. The new 
feature is the presence of two new poles located at p ^ m£)/2 ± m^ i corresponding to a tachyon-like state with mass 
~ ms = Mp/y/a. It is easy to prove that this tachyon is a ghost. The residue of the pole is complex and thus the 
norm of the state is not positive. But this ghost is very different from that arising in the self-accelerating branch of 
the usual DGP model. In the latter case the ghost couples gravitationally to brane-localized matter and has mass 
~ niD ~ Hq^ so it can be excited in the IR, exactly at the scale where we expect significant modifications from the 
self-accelerating branch solution. 

Conversely in the DGP model with higher curvature brane terms the ghost is very heavy, it has mass ~ Mp / ^/a. 
Despite of its large mass, the tiny real part of tachyon pole leads to an instability of the Minkowski spacetime with 
time-scale of order {mu)~^- Self- accelerating solutions exist even in the ordinary branch of the DGP model if one 
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FIG. 1: Schematic view of the Euclidean pole structure for the DGP model with higher curvature brane terms in the ordinary 
branch and for q > 0. The massless particle does not propagate any degree of freedom. This is a consequence of the absence 
of normalizable massless 4D states in the model. The pole on the real axis corresponds to a resonance with a typical mass 
~ ran ~ lO~'^^GeV . The resonance does not lead to any instability of the Minkowski spacetime. The two complex poles 
correspond to a massive tachyon, m ~ Mp/y^ ~ W^^GeV. 



considers higher curvature terms on the brane. Furthermore, since in the ordinary branch the gravitational interaction 
is mediated by the resonance of mass mo, we naturally expect that at large distances the heavy particle decouples 
from brane-localized matter. 

The physical explanation of this behavior involves both 4D brane and 5D bulk physics. The tiny real (in the 
Euclidean space) contribution, of order mo, to the pole describing the ghost has to be explained in terms of 5D non 
local effects, hence as a IR effect. On the other hand the huge imaginary contribution of order Mpj ^fa has to be 
explained has a short distance 4D effects which modifies only the UV behavior of the propagator. The net result is 
a very massive tachyon which has a very small and positive real part that, as in the self-accelerating branch of the 
usual DGP model, leads to an instability of the Minkowski space. This suggests that in the DGP model with higher 
curvature brane terms even the ordinary branch is self-accelerating but without some of the pathologies arising in the 
self-accelerating branch. 

The essential features of the model discussed here seem quite general and in principle are not limited to brane- 
world scenarios. They could be present, at linearized level, in any theory of gravity whenever the pole structure 
of the Euclidean propagator of some excitation is determined by both IR (positive real part) and UV (imaginary 
part) effects. The two mass regimes are so disentangled that one can have a heavy excitation with mass of order 
A/Re(p)2 -I- Im(p)2 ^ Re(p) leading to a tiny instability of the Minkowski background. 

The main difficulty that we have to face in order to give a consistent physical explanation of the above-described 
instability of Minkowski vacuum is the ghost nature of the excitation. The corresponding pole has a complex residue 
and this leads to a non-unitary theory. Obviously this unitarity problem cannot be solved in the framework of a low- 
energy effective theory of gravity described by the action One can naturally assume the UV completion of our 
theory of gravity to be ghost free. The ghost could be reabsorbed in the (infinite) series of terms of the fundamental 
theory. Whether the instability due to the ghost does or does not remain in the fundamental theory is a subtle issue 
that it is not addressed here. 

The presence of higher curvature term on the brane seems to lead quite generically to an instability of the Minkowski 
spacetime in the ordinary branch. There are indications that adding a /(i?)-term on the brane leads to self- accelerating 
solutions even in the ordinary branch [26| . Interestingly enough, there is no trace of this cosmological instability at 
the linearized level because the propagator remains unchanged by adding /(i?) brane terms. In this case instabilities 
of the Minkowski spacetime arise at the full non-linear level only. This is another example in which nonlinear effects 
in the DGP model have a strong impact on the dynamics. The study of cosmological solutions arising from action 
([5]) will be presented elsewhere. 

Although it is not a central issue, let us we briefly discuss the case — ao < a < 0. Now the solutions of Eq. ((TU]) (with 
the -f sign) are 



X\ = mo 



1 



-arctan 



1 



y 



>o, 



a;2,3 



mo 
2 



tan ( —arctan 



y 



<0, (16) 
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where y = —a/ao < 1. In this region we find two resonances and a tachyon, all of them with zero imaginary part. The 
pole structure is richer than in the usual DGP case, but no new features arise in the ordinary branch. The behavior 
in the self-accelerating branch is qualitatively the same as that described above. In particular, the poles structure 
in the Euclidean plane is y-reflected. This means that, going from one branch to the other, we just need to change 
resonances into tachyons and vice versa. This general result holds in the usual DGP case too. Always a very light 
tachyon- like state appears in this branch and it is a ghost with a mass of the order mo- 

Higher curvature brane terms do not change the large distance behavior of the static potential of the usual DGP 
model. Nevertheless they can give short distance contributions, at energy scales of order Mp/y/a. The modified DGP 
model ^ corrects GR both in the IR and in the UV. Neglecting the tensorial structure of the propagator ([9]), the 
static potential is 



Vir) 



dt 



(17) 



where 4>{p, ?; = 0) is the Fourier-transformed scalar part of the propagator ^ on the brane. Focusing on the 
phenomenologically interesting case a > the potential reads 



Vir) 



1 



(2Vi(r)- Re[(2 + 3i77)T/2(r)]) , 



where 77 — mi/m2 and Vi, V2 are obtained from the standard DGP potential [ij 



VDGp{r) = 



sin(m£)r)Ci(m£)r) + — cos(m£)r)(7r — 2Si(m£)r)) 



(18) 



(19) 



substituting 
becomes 



mi and mo ^ mi /2 + 11121 respectively. In the large-r limit (r ^ 1/mi, l/m2) the potential (fT5)) 



V{r) 



1 



miaiT^{4 + 977^) 



2-37? 

{m2rY 



m2r 



2 cos 



mir 



Srj sin 



mir 



(20) 



As expected, corrections due to higher curvature brane terms decay exponentially. As in the usual DGP model 
the gravitational potential is five dimensional and its behavior can set a limit on the DGP parameter, mjj. In the 
intermediate region, 1/7712 <C t' ^ 1/771.1 the potential is 



1 / 1 277 , , . ^^ 377 1 

77^:^ + ^ - log(mir)) + ^- ^ 

77i2Q!7r(4 + 97;"^) \m2r ir tt [m2ry 



2e- 



3772 



7712 r 



(21) 



where 7 is the Euler constant. In this region one recovers Newton potential corrected by subleading logarithmic terms 
arising from the DGP bulk and by exponential and by corrections arising from the higher curvature brane terms. 
In the small-r region, r <C 1/7772, l/^i, the potential behaves as ~ r. This, together with submillimeter table-desk 
experiments can put constraints on the value on the parameter a. 

In this paper we have investigated the DGP model corrected by higher curvature brane terms. We have found that 
these corrections have a dramatic effect on the linearized spectrum of the model. Ghost excitations with mass of 
order Mp, which describe an instability of Minkowski vacuum with time-scales of order H^^^aie generated. In the IR 
the tachyonic excitation is expected to decouple from brane-localizcd matter. Our modified DGP model represents 
therefore a very promising framework for solving the cosmological constant problem, in which Planck scale physics is 
responsible for the elementary excitations driving the accelerated expansion of the universe, but the timescale of the 
instability is settled by the gravitational physics at large scales. 

Our result gives a strong hint that short distance corrections to the DGP model may held the key for solving some 
of the open problem of the model. There are still some issues and consistency checks, which have been not addressed 
in detail in this paper. The first is the particle physics meaning of the ghost excitation driving the instability of 
Minkowski space. The main issue here is not the ghost nature of the excitation - it is rather natural to assume 
that this problem will be solved by the UV completion of the theory, because any consistent quantum theory of 
gravity, such as string theory, is expected to be unitary- but the reliability of our low-energy approximation given 
by the action ([5]). It is not clear to what extent the truncation is consistent and if in an alternative treatment the 
ghost excitation would still be present. An other, related, issue is the reliability of perturbation theory. Our results 
have been derived in the linear approximation. It is well known that the linear approximation of the DGP model is 
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strongly limited because of the vDVZ discontinuity and of the strong coupling effect. Nonlinear effects are therefore 
crucial for making the model consistent. It is not clear whether the features we have found at the linearized level 
still persist at a full nonlinear level. A very important check in this context is represented by the investigation of 
the cosmological solutions of the our model. It is also important to stress that the introduction of higher-derivatives 
interactions -also in the form of higher curvature terms- in the DGP model may be also very useful for solving the 
vDVZ discontinuity problem. These terms can support a robust implementation of the Vainshtein effect [l^ and still 
allow for self-accelerating de Sitter solutions with no light ghost instabilities [2^. 
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